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Abstract. This paper is the continuation of a previous one [L. Samaj and B. 
Jancovici, 2007 J. Stat. Mech. P02002]; for a nearly classical quantum fluid in a 
half-space bounded by a plain plane hard wall (no image forces), we had generalized 
the Wigncr-Kirkwood expansion of the equilibrium statistical quantities in powers of 
Planck's constant h. As a model system for a more detailed study, we consider the 
quantum two-dimensional one-component plasma: a system of charged particles of one 
species, interacting through the logarithmic Coulomb potential in two dimensions, in 
a uniformly charged background of opposite sign, such that the total charge vanishes. 
The corresponding classical system is exactly solvable in a variety of geometries, 
including the present one of a half-plane, when (3e 2 = 2, where (3 is the inverse 
temperature and e is the charge of a particle: all the classical n-body densities are 
known. In the present paper, we have calculated the expansions of the quantum 
density profile and truncated two-body density up to order fr 2 (instead of only to order 
h in the previous paper) . These expansions involve the classical n-body densities up to 
n = 4, thus we obtain exact expressions for these quantum expansions in this special 
case. 

For the quantum one-component plasma, two sum rules involving the truncated 
two-body density (and, for one of them, the density profile) have been derived, a long 
time ago, by heuristic macroscopic arguments: one sum rule is about the asymptotic 
form along the wall of the truncated two-body density, the other one is about the 
dipole moment of the structure factor. In the two-dimensional case at /3e 2 = 2, we 
have now explicit expressions up to order h 2 of these two quantum densities, thus 
we can microscopically check the sum rules at this order. The checks are positive, 
reinforcing the idea that the sum rules are correct. 
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1. Introduction 

The model under consideration is the two-dimensional (2D) one-component plasma (also 
called jellium). This model consists of one species of charged particles in a plane. Each 
particle has a charge e and a mass m. Two particles, at a distance r from each other, 
interact through the 2D Coulomb interaction v(r). This interaction is determined by the 
2D Poisson equation V 2 w(r) = — 2ne 2 S(r), the solution of which is v(r) = — e 2 ln(r/ro), 
where ro is an arbitrary length which only fixes the zero of this potential. In addition, 
there is a charged uniform background of charge density opposite to the particle charge, 
so that the total system is neutral. 

When the inverse temperature f3 is such that the dimensionless coupling constant 
(3e 2 = 2, the equilibrium statistical mechanics of the classical (i.e. non-quantum) system 
is completely solvable in a variety of geometries, in particular when the background and 
particles are confined into a half-space by an impenetrable rectilinear plain hard wall 
(there are no image forces): all the classical n-body densities are known [1, 2J. 

In its three-dimensional version, the quantum one-component plasma is not only 
of academic interest. It has been used, in first approximation, as a model for the 
electrons of a metal |3]. In the present geometry of a half-space, this model might be 
used for describing what happens near the surface of the metal. Since the logarithmic 
interaction v{r) is the Coulomb potential in two dimensions, the 2D one-component 
plasma is expected to have general features which mimic those of the 3D one. 

Under certain conditions, the equilibrium properties of an infinite quantum fluid, 
in the nearly classical regime, can be expanded in powers of Planck's constant h: this 
is the Wigner-Kirkwood expansion [HE]. In a previous paper [6], we have generalized 
the Wigner-Kirkwood expansion to the case of a quantum fluid occupying a half space. 
We have obtained expressions for the first quantum correction of order h to the density 
profile and to the two-body density, in terms of some n-body densities of the classical 
fluid. In section 3 of the present paper we extend these calculations to order h 2 . Actually, 
instead of h, we use the thermal de Broglie wavelength proportional to it, A = hJ p/m. 



For the quantum one-component plasma in a half space, by heuristic macroscopic 
methods, some sum rules involving the one-body density profile and the two-body 
density have been derived. The purpose of the present paper is to check these sum 
rules at order h 2 , in the special case of the 2D one-component plasma at (3e 2 = 2, using 
the generalized Wigner-Kirkwood expansion, which is a microscopic approach. 

The paper is organized as follows. Section 2 brings a recapitulation of the method 
for constructing the expansion of the quantum Boltzmann density in configuration 
r-space for fluids constrained to a half-space [6] . This expansion is subsequently used in 
section 3 to compute the quantum one-body profile and two-body density to order h 2 
for the studied 2D one-component plasma. In section 4, three sum rules are reviewed: 
a perfect screening rule, a sum rule about the asymptotic form of the two-body density 
along the wall, and a dipole sum rule; their expansions to order h 2 are given. In section 
5, some classical n-body correlation functions, which are needed, are studied. In section 
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6, the perfect screening sum rule is shown to hold, at order h 2 , for the quantum 2D one- 
component plasma at f3 — 2. In sections 7 and 8, the same is done for the asymptotic 
form and the dipole sum rules, respectively. Section 9 is a Conclusion. 



2. Boltzmann density for the half-space geometry 

We first consider a general quantum system of N identical particles j — 1, 2, . . . , N of 
mass m, formulated in v space dimensions. Particle position vectors r 1; r 2 , . . . , are 
confined to the half-space A defined by Cartesian coordinates r = (x > 0, r- 1 ), where 
r 1 - G R" _1 denotes the set of {y — 1) unbounded coordinates normal to x. As usual, 
we start with a finite A" and a finite volume |A|, and later we take the thermodynamic 
limit A" and |A| going to infinity (when A becomes a half space) with a finite mean 
number density n = N/\A\. The hard wall in the complementary half-space A of points 
r = (x < 0, r -1 ) is considered to be impenetrable to particles, i.e. the wavefunctions 
of the particle system vanish as soon as one of the particles lies at the wall. For the 
sake of brevity, we denote the i/A^-dimensional position vector in configuration space by 
r = (r 1; r 2 , . . . , r N ) and the corresponding gradient by V = (Vi, V 2 , . . . , Vjv)- In the 
absence of a magnetic field, the Hamiltonian of the particle system is given by 

H=±-(-ihv) 2 + V{r), (2.1) 

where h stands for Planck's constant and ^(r) is the total interaction potential. 

For infinite (bulk) quantum fluids of particles interacting via pairwise sufficiently 
smooth interactions with neglected fermion/boson exchange effects, Wigner [4] and 
Kirkwood [5j constructed a semiclassical expansion of the Boltzmann density in 
configuration space (at inverse temperature f5), -Bg(r) = (r|e _/3H |r), in even powers 
of the thermal de Broglie wavelength A = h{(5 /m) 1 / 2 . Recently [6], we have generalized 
the Wigner-Kirkwood method to quantum fluids constrained to the above defined half- 
space A. The final result for the Boltzmann density in configuration space was obtained 
as a series 

oo 

B,{v) = Y J Bf\v), (2.2) 

n=0 

where the terms B^\v) with n = 0, 1, 2, . . . can be calculated systematically with the 
aid of an operator technique. 

The result for Bp (r) was found in the form 

i N 

B ^ = ym^^ v ^- e ' 2xV ^- (2 - 3) 

The bulk counterpart of this term corresponds to the classical Boltzmann density 
oc e _/3y . Here, each particle gets an additional "boundary" factor 1 — exp(— 2x 2 /X 2 ) 
which goes from at the x = boundary to 1 in the bulk interior x — > oo on the 
length scale ~ A. The product of boundary factors then ensures that the quantum 
Boltzmann density vanishes as soon as one of the particles lies on the boundary. The 
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dependence of the boundary factor on the de Broglie wavelength A is non-analytic; this 
fact prevents one from a simple classification of contributions to the Boltzmann density 
according to integer powers of A like it is in the bulk case. However, when in the 
calculation of statistical averages the exponential part exp(— 2x 2 /A 2 ) of the boundary 
factor is integrated over the x-coordinate, the analyticity of the result in the parameter 
A is restored. At this stage we only notice that when the product of boundary factors 
is expanded as follows 

N N i N 



n(i- 

i=i 



-2x 2 /X 2 



i-E< 



-2x 2 /X' 



+ ±Y: e- 2 ^ A V 2 ^ A2 + • • ■ , (2.4) 



j,k=l 



the integration of each exponential term exp(— 2x 2 /X 2 ) over x produces one A-factor as 
the result of the substitution of variables x = Xx'. 



The result for B 



(i), 
P 1 



r) reads 
1 



N N 



^X) 



E n (l-e- 2 ^ /A > fc e- 



k=l 3=1 



d_ 

dxi 



-pv 



+ e 



N 

-«< JJ 

3=1 



-2x 2 /X' 



A 



6 



;2.5) 



Here, the dependence on A appears also via the combination x exp(—2x 2 / X 2 ) . This 
function has a maximum of order A and therefore it is a legitimate expansion parameter. 
When integrated over the particle coordinate x, it gives a contribution of order A 2 , 
"weaker" than A. 

Keeping all contributions up to the order A 2 in the Boltzmann term ai 2 (r), one 



has 



1 



N 

n(i 



-2x 2 /\' 



Note that, with regard to the equality 



A 



D 



+ o(A 2 



V e 



-pv 



-pv 



(3 2 (VV) 2 -PV 2 V 



;2.6) 



(2.7) 



one can eliminate the squared gradient term in favour of the Laplacian term in equations 
(1231) and (EE}. 

The first three Boltzmann terms (12.31) . (12.51) and (I2.6P exhibit properties analogous 
to their bulk counterparts: the maximum of B^(r) is of order X n . We anticipate 
that this formal structure is also maintained on higher levels. As a consequence, the 



knowledge of the first three Boltzmann terms Bg, Bg and B Z> is sufficient in order 



? (2) 



11 



to obtain the expansion of the quantum Boltzmann density up to the A 2 order: 



^(r) = Bf{v) + Bf{f) + BV{1) + o(A 



?(2)/ 



'P K*J~r*Jp ^j-TJ^g y±j-ru V ij. (2.8) 

The quantum fluid of present interest is the one-component plasma (jellium) in 
v = 2 space dimensions. The system is composed of iV mobile pointlike charges e, 
neutralized by a uniform oppositely charged fixed background. The total interaction 



Correlations and sum rules in a half-space quantum plasma 



5 



potential V(f) satisfies for each of the particle coordinates the Poisson differential 
equation 

N 

V 2 V(v) = -2ne 2 £ Sfa - r k ) + 2ne 2 n, j = 1, 2, . . . , N. (2.9) 



fc=i 



Here, the second term on the right-hand side (rhs) comes from the particle-background 
interaction and n = iV/|A| is the mean number density of the mobile charges. The 
summation over the particle index j of the set of N Poisson equations (12.91) results in 

^2 N 

V V(r) = -27re 2 ^ Sfa - r k ) + 2ne 2 Nn. (2.10) 



3,fc=l 



The first term on the rhs of (I2.10p . when weighted by the classical Boltzmann 
factor oc Y\j<k \ r j ~ ^kl 13 ^ which vanishes at zero interparticle distance, does not give 
any contribution to the Boltzmann density. The application of the equality (12.71) and the 
replacement of V V(r) by the particle-background term 2nNe 2 n in the relations (I2.5D 
and (12. 6p simplifies substantially the expansion formula for the quantum Boltzmann 
density fl2T8|) : 

Bp (f) = , J , J e~ pv ( 1 - ^ 2K0e 2 Nn) T7 (l - e~ 2x V^' 



-± L-PV ( i _ ^-2n(3e 2 Nn) TT (l - e 



N N q 



+ E II (l-e~ 2x V^) Xk e- 2 ^ 2 ^e^ v 



+ n(l-e-^ /A2 )^Ve-^| + o(A 2 ). (2.11) 
3. Statistical quantities for the half-space geometry 

According to the standard formalism of statistical quantum mechanics, the partition 
function of the iV-particle fluid (with ignored exchange effects) is given by the integration 
of the quantum Boltzmann density over configuration space: 

The quantum average of a function /(f) is defined as follows 

</>au = ^4 / df ^(f)/(f). (3.2) 



At the one-particle level, one introduces the particle number density 

M r ) = (f>( r - r ;)) ■ ( 3 - 3 ) 

\ j=l I qu 

At the two-particle level, one considers the two-body density 

< 2 u V,r') = ( E 5(r-r,)5(r'-r k )) . (3.4) 

\ 3,k=l I qu 



Correlations and sum rules in a half-space quantum plasma 



6 



It will be useful to consider also the truncated two-body density 

™qu )T Oi> r 2 ) = 7ig)(ri, r 2 ) - n qu ( ri )n qu (r 2 ) (3.5) 

vanishing at asymptotically large distances |ri — r 2 | — > oo. The general multiparticle 
densities are defined in analogy with (13.41) . i.e. the corresponding product of 5-functions 
is summed out over all possible multiplets of different particles. 

The classical partition function Z and the classical average of a function /(r) are 
defined as follows 

^-iw.Lyixr^^ (3 ' 7) 

The classical values of statistical quantities will be written without a subscript, like 
n(r), (r, r'), etc. In the calculations which follow, we shall need explicitly truncated 
forms of the classical three-body density 

n (3)T (ri, r 2 , r 3 ) = n (3) (ri, r 2 , r 3 ) - n (2)T (ri, r 2 )n(r 3 ) - n (2)T (ri, r 3 )n(r 2 ) 

- n (2)T (r 2 , r 3 )n(r x ) - 7i(r 1 )n(r 2 )n(r 3 ) (3.8) 

and of the classical four-body density 

n (4)T (ri, r 2 , r 3 , r 4 ) = n (4) (n, r 2 , r 3 , r 4 ) - n (3)T (ri, r 2 , r 3 )n(r 4 ) 

- n (3)T (ri, r 2 , r 4 )n(r 3 ) - n (3)T (ri, r 3 , r 4 )n(r 2 ) 

- n (3)T (r 2 , r 3 , r 4 )n(ri) - n (2)T (ri, r 2 )n (2)T (r 3 , r 4 ) 

- n^(r u r 3 )n( 2 ) T (r 2 , r 4 ) - n( 2 ) T (r 1; r 4 )n( 2 ) T (r 2 , r 3 ) 

- ra (2)T (ri, r 2 )n(r 3 )n(r 4 ) - n ( - 2)T (r 1 , r 3 )n(r 2 )n(r 4 ) 

- n (2)T (ri, r 4 )n(r 2 )n(r 3 ) - n (2)T (r 2 , r 3 )n(r 1 )n(r 4 ) 

- n (2)T (r 2 , r 4 )n(ri)n(r 3 ) - n (2)T (r 3 , r 4 )n(ri)n(r 2 ) 

- n(ri)n(r 2 )n(r 3 )n(r 4 ). (3.9) 

In what follows, we shall restrict ourselves to the model system of the one- 
component plasma constrained to the two-dimensional half-space A. Now, the Cartesian 
coordinates of r become (x,y), with the origin on the rectilinear plain hard wall, the y 
axis along the wall, and the system occupying the x > half-space A. 

3.1. Partition function 

Substituting the expansion of the Boltzmann density (12.111) into the definition (13. ip of 
the quantum partition function Z qu and performing expansions of type (12 .4p for the 
products of boundary factors, we obtain 



Z qu _ ^ A 2 

~Y " 24 

+ ' 



27c(3e 2 Nn- f dne^^fifn) 

J A 

\ / drx / dr 2 e- 2 ^/ A2 e- 2 ^/ A2 n( 2 )( ri ,r 2 ) 

2! J A J A 
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J A 



d 



X 2 



d 2 



^^ n ^ + 24A dri ^! n(ri)+ ° (A2) - (3 - 10) 

Here, we keep in mind that the integration of an exponential term exp(— 2x 2 /A 2 ) over 
x produces one A-factor. 



3.2. One-body density 



To calculate the quantum one-body density (13.31) . we take advantage of the invariance 
of the Boltzmann density (12. lip with respect to permutations of the particle indices and 
write down 

' ' 1 N rJ A & B f>W8(T-*i)- (3-11) 



n qu {r) 



Z qu N\ 



In each term of the Boltzmann density, we separate the "reference" ^-dependent part, 
which is kept unchanged, and expand in analogy with (12 Ah the remaining part dependent 
on (r 2 , . . . , r N ) coordinates. Like for instance, 

N ( N 



1-Ee 

i=2 



-2x 2 /\ 2 



N 



+ 



2! 



e -2xj/\ 2 e -2xl/\ 2 



(3.12) 



j,k=2 



The quantum partition function Z qu in the denominator on the rhs of ( 13. lip is 
substituted by the expansion ( I3.10P and subsequently expanded in "virtual" A powers. 
After simple but lengthy algebra, one obtains 



Wqu(ri) 



n r 



1 - e" 2 ^ 

H — - / dr 2 / dr^ e 

2! J a 2 Ja 3 



dr 2 e- 2 ^ A2 n^ T ' 



ri,r 2 ) 



^i,r 2 ,r 3 ) 



+ / dr 2 e" 
J A 



2x\l* 



X 2 



d 

dx 2 



+ 



A 2 



24 



U 2 



dx 2 



d 2 



dl "2 TT2 n 

dx{ 



(2)T, 



ri,r 2 ) 



+ e 



-2x\l* 



d 



n{r l 



f d^e-^/^n^^,^) 

J A 



(3.13) 



Let us assume that the classical averages under integrations in (13.131) are analytic 
functions of the x-coordinate at the boundary x = 0. For instance, the classical 
density profile n(r) = n(x) is supposed to exhibit the Taylor expansion n(x) = 
n(0) + n'(0)x + n"(0)x 2 /2\ H . Then, the integral 



dxe- 2x2/x2 n(x) = A / dx' e~ 2x ' 2 n(\x') 
o Jo 



A^M0) + A 2 ^'(0) + O(A 3 ). 



(3.14) 
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Performing an analogous procedure in all integrals on the rhs of (13.131) . we finally arrive 
at the expansion of the one-body density up to the A 2 order: 

rv( ri ) = (l - e" 2 ^ 2 ) |n(n) - A^|/ d^n^n, (0,y 2 )] 
+ y (|) / d ^ / dy 3 n^ T [r 1; (0,y 2 ), (0,^)] 



A 2 
+ 24 



«9 2 /■ 

^! n(ri) - y ^ 



dn( 2 > T [ri, {x 2 ,y 2 )) 



x 2 =0. 



+ e -2^/A 2 



d 



dx 2 

X^jdy^^^y,)} 



(3.15) 



3.3. Two-body density 



To calculate the A-expansion of the quantum truncated two-body density given by 
relations H3.4[) and (13.51) . we proceed as in the previous subsection. We first take 
advantage of the permutation invariance of the Boltzmann density (12.111) to write down 

1 N(N-l) 



^qu 



dlBp(f)6(T-T 1 )6(T f -T 2 ). 



(3.16) 



Then we separate the reference part in the Boltzmann density (12.111) which depends on 
coordinates rj., r 2 and expand in analogy with ( 13.121) the remaining part. Finally, using 
the integration procedure of type (13. 14f) in all integrals we arrive at the expansion of 
the two-body density up to the A 2 order: 



(2)T 
qu 



(r 1; r 2 )= (l-e^ 2 ) (l - e" 



■2x 2 /X 2 



n< 2 > T (n,r 2 ) 



-X^jdy 3 n^[v u r 2 ,(0,y 3 )] 

dy 3 J dy A n^ T [T U T 2 ,(0,y 3 ), (0, y 4 )] 



+ 



A_ 2 /7T 

~2 

A 2 



24 

+ 



d 2 d 2 
+ 



^ dx 2 dx\ 

^-2x 2 /X 2 



n (2)T (n,r 2 ) 



dy; 



0ra (3)T [ri,r 2 , (x 3 ,y 3 )} 



dx 3 



x 3 =0- 



2 l\2 



xie 



-2x\/X 



dx\ 



n( 2 > T (r 1; r 2 ) 



-X^jdy, ^[^,^,(0,2/3)]} 



-2x1/ X 2 



x 2 e 



-2x 2 /X 2 







- A 



dx 2 y 

J dy 3 n^ T [r u r 2 ,(0,y 3 )] 



n (2)T (r 1; r 2 ) 



(3.17) 
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4. Review of the sum rules 

The following sum rules [except the perfect-screening one (14.11) ] apply only to the one- 
component plasma (not to many-component ones); they rely on the facts that, for 
the one-component plasma, the mass and charge fluctuations are proportional to one 
another and the static resistivity vanishes. Although the sum rules can be written for 
a quantum one-component plasma in a i/-dimensional half-space with any v, here we 
consider only the case v = 2. Some of the sum rules that we review were originally 
written more generally for the time-displaced correlations, but here we consider only 
their static limits. 

The perfect screening sum rule expresses that the charge cloud around a particle of 
the system has a charge opposite to the charge of this particle. This sum rule has the 
same form as in the classical case: 

dri7^ )T (ri,r 2 ) = -n qu (x 2 ). (4.1) 

Although we have no doubt about the validity of (14. ip . for a check of the calculations in 
section 3 and of their application to the present 2D one-component plasma at (3e 2 = 2, 
we shall verify ( 14. ip at order A 2 in section 6. 

Macroscopic arguments gave the asymptotic form of the quantum truncated two- 
body density along the wall [7j 

(2)T/ \ / \ X li X V I a q\ 

<J ( r i, r 2), ~ 7 rj, (4.2) 

lm-M|-»°o [yx - y 2 ) 1 

(after perhaps an averaging on local oscillations in yi — y 2 ) with the sum rule 

oo roo | 

dxi / dx 2 f(xi, x 2 ) = — 5 \2tiu) B coth((3tiu) B /2) — huj p coth.((3hu p /2)] (4.3) 
o Jo ^K l e A 

(here j3e 2 = 2), where the bulk and surface plasma frequencies, for two dimensions and 

a plain hard wall, are, respectively, 

\ m J \ m J 

Any microscopic check of (14.31) would be welcome. The expansion of the r.h.s. of (14.31) 
in powers of H starts with the classical value — l/(47r 2 ) and the next term is of order 
h 4 . Therefore, at order A 2 , in section 7 we shall only be able to check that there is no 
quantum correction. 

Finally, another macroscopic argument gave the quantum form of the dipole sum 
rule [8] 



roo r poo r 

dx 2 dxixi J dyi^ 2)T (ri,r 2 ) + x 2 n qu (x 2 ) 



hu p (3fao- 



,coth- — v - (4.5) 
4vre 2 2 v ; 

(here j3e 2 = 2). The quantity between square brackets in (14.51) is the dipole moment of 

the structure factor. Here too, a check would be welcome. Now, the expansion of the 

rhs of ( 14.51) in powers of h starts with the classical value — l/(47r) and the next term is 

— A 2 /(127ra 2 ), where a is the average interparticle distance defined by n = l/(ira 2 ). In 

section 8 we can check this sum rule at order A 2 . 
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We shall need some information about the classical densities of the 2D one-component 
plasma at f3e 2 = 2. We express all lengths in units of the average interparticle distance 
a, thus n = 1/tt. The density profile is p] 

2 



n(x) = n 



IT JO 



°° df exp[-(* - xV2) 2 ] 



1 + $(t) 



(5.1) 



where $(t) is the probability- integral function. The two-body truncated density is 



n 



(2)T. 



ri, r 2 ) = -n 2 exp[-2(ar 2 + a^M 1 "!) r 2)k(r 2 , n), 



where 



fe(ri,r 2 ) 



oo exp 
dt 



-t 2 + t{x! +x 2 )V2- it{y x - y 2 )V2 



(5.2) 



(5.3) 



TWO l + $(t) 

Higher-order ra-body truncated densities contain a sum of product of n factors k, the sum 
running on all oriented cycles built with i"i,r 2 , ...,r n [2]. One finds for the three-body 
truncated density 



71 



(3)T. 



ri, r 2 , r 3 ) = n 3 exp[-2(x 2 + x 2 + x 2 )] [fc(ri, r 2 )fc(r 2 , r 3 )/c(r 3 , r x ) 



+ /c(ri,r 3 )A;(r3,r 2 )A;(r 2 ,r 1 ) 
and for the four-body truncated density 



(5.4) 



71 



(4)T, 



ri, r 2 , r 3 , r 4 ) = - n 4 exp[-2(a;? + x\ + x 3 + x 2 )] [fc(ri, r 2 )£;(r 2 , r 3 )/c(r 3 , r 4 )fc(r 4 , r a ) 
+ fc(r!, r 2 )fc(r 2 , r 4 )fc(r 4 , r 3 )fc(r 3 , n) 
+ fe(ri, r 3 )A;(r 3 , r 2 )£;(r 2 , r 4 )fc(r 4 , n)] 

+ complex conjugate. (5.5) 

We shall also need some integrals of these densities. In f)5.2p . the product of the k 
functions contains J °° dt exp(— ityi\f2) J °° dt' exp (it'?/! a/2), therefore the integral of this 
quantity on yi is ix\[2 5(t — t'). One obtains 

2t 2 + t(xi + x 2 )2 v / 2 



dyin^'^ri, r 2 ) = -// 

By the same method, one finds 

3 8^2 



roo exp 

4V2exp[-2(:r 2 + x 2 2 )} / dt 



[1 + «&(*)]* 



(5.6) 



dy 3 n 



(3)T. 



ri,r 2 ,r 3 ] 



n 



x 



x 



7T 



exp[-2(x 2 + x 2 + x 2 )] 



°° ^ exp[-t 2 + t(x! + x 2 )y / 2 - it(yi - y 2 ) ^2] 
1 + $(t) 

°° d ^ / exp[-2t /2 + t / (x 1 + x 2 + 2x 3 )V2 - it'(y 2 - Vl )V2] 



[1 + *(tf)] s 



+ complex conjugate, 



dyi J dy 3 n (3)T [ 



ri,r 2 ,r 3 ] 



n 



! 320f exp[-2(x 2 + x\ + x^)] 



dt 



exp[-3t 2 + t(xi + x 2 + o: 3 )2a/2] 
[l + $(t)]3 



(5.7) 



(5.8) 
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dy 3 J d^ (4)T [r 1 ,r 2 ,(0,|/3),(0,i/4)] = -n 4 32exp[-2(^ + ^)] 

roo 
X / dt' 

Jo 

+ / dt 
Jo 



dt i + m 






,exp[-3t' 2 + t'ixi +x 2 )V2~ 


- it' (2/2 


-yi)V2] 


[1 + $(t')] 3 






exp[-2t 2 + t(x l + x 2 )y/2-: 


Hvi - 


V2)V2] 


[l + $(t)P 






^,exp[-2t' 2 +t'(x 1 +x 2 )V2- 


- i*'(2/2 


-yi)V2] 







[i + m? ^' " / 



+complex conjugate, (5.9) 

| dj/i | dy 3 / d^n^fn,^, (0,1/3), (0,y 4 )] = -n 4 1927rv / 2exp[-2( a : 2 + a ; 2 )] 
/°° exp[-4t 2 + t( 3 ; 1 +x 2 )2v^] 

x /o dt [TTw • (5 - 10) 

6. Perfect screening 

Omitting some terms which do not contribute to the sum rules at order A 2 , and using 
the results in section 5, we find from (3.17) 

qu v n,r 2 ) = [1 -exp(-2:r 2 /A 2 )][l - exp(-2x 2 2 / X 2 )} 

exp[-3t 2 + t(x! + x 2 )2%/2\ 



-Xn^fdt* [i + m]3 

Jo [1 + *(t)] 4 



AVl2V2vr 2 /-dt- eXp [- 4t2 + ^ +a;2)2 ^ 
— exp[-3t 2 + t(xi + 2:2)2^] 1 



-X 2 n 3 -y/2% I dtt 



3 v ^y„ [i + « W ]' 7 

+ (" A2 ^ (S + - n 2 4V2e X p(-2,;/A 2 ), 1 9 



6 \(9x 2 (9^2 / ~ 9xi 

-) 

dx 2 J 

00 exp[-2t 2 + t(x! + ^2)2^] 1 



-n 2 4V2exp(-2a%/\ 2 )x 2 - 



On the other hand, we find from (3.15) 

n 9u (x 2 ) = [1 - exp(-2^/A 2 )] | A ne xp(-2x^) J~ ^'f+^f^ 



2 1 d 2 

+ \2^dx\ 
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, ,„ , 9x f°° , exp(-3t 2 + tx 2 2\/2) 

+ xw"eM-24)l ^ pl |1 + t(t)] 2 3 

exp(-2^) fdt^^ 

FV 2; 7o i + $(t) 

+An-exp(-2x 2 )y o dff j 

, „ , m9n d { 2 , 9x Z" 00 , exp(-t 2 + te 2 2\/2) 
+ exp(-2^/AV 2 - exp(-2,i) j£ dt ; - ^ 

+ A^2^exp(-24) ^ dt exp(-2^^2V2) | + . . . (M) 

Using [9] 

/ dxexp(-2a; 2 + ta;2v / 2) = W-exp(t 2 )[l + $(t)], (6.3) 
jo V 8 

we can compute the integral on x\ of ( 16. lj) at order A 2 and check that it is equal to the 

opposite of (16.21) . 

The sum rule (14.11) is verified. 

7. The sum rule about an asymptotic form 

For obtaining the asymptotic forms of the classical rz-body densities as \y\ — 2/2 1 — > 00, 
one uses the integration per partes for the Fourier transform of a function F(t) 

rdtF(t)exp[-ityV2]~^. (7.1) 

The different classical n-body truncated densities or their integrals which appear in 
(3.17) are all found to have an asymptotic form oc exp[— 2(x 2 + x\)}/{yi — y 2 ) 2 - This 
gives the asymptotic form (14.21) . where (some terms which do not contribute to the sum 
rule at order A 2 have not been kept) 

f(x u x 2 ) = [1 - exp(-2a; 2 /A 2 )][l - exp(-2a: 2 /A 2 )] 

x exp[-2(x 2 + x 2 )] (~™ 2 | - An 3 4 - AV6tt) + • • • . (7.2) 

Now, we note from equation (13.141) that, at order A 2 , an integral of the form 
/ oo dxexp(-2a; 2 /A 2 )F(a;) becomes A v /7r/8F(0) + A 2 (l/4)F , (0). Then, from (Q, taking 
into account that in our units n = l/V, 

/•oo roo 1 / 1 1 1 \ 

Jo dXl Jo dX2fiXl > X2)= -^ + X {^ + ^-2^) 

+ A2 (-i + i + i-i) + ° (A2) 
' «(A 2 ) (7.3) 



4tt 2 

where A is in units of a. 

At order A 2 , there are no quantum corrections, in agreement with section 4. 
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8. Dipole sum rule 

Omitting some terms which do not contribute to the sum rule at order A 2 , using [9] 

J dxxexp(-2x 2 + 2v / 2te) = - + ^texp(t 2 )[l + $(t)], (8.1) 

and noting that, at order A 2 , an integral of the form / °° dxx exp(— 2x 2 / X 2 )F(x) becomes 
(1/4)A 2 F(0), one finds from (16. ip . contributing to the sum rule at order A 2 , 

[1 - exp(-2x 2 2 /X 2 )] exp(-2x 2 2 ) 

-n 2 \[2 



dxiXi / dyin^ 2)T (ri,r 2 ) 



dt 



[1 + $(t)] 2 

j_\2 2 /o /■°°^ exp(-2t 2 + tx 2 2V2) 

+A n V2 / dt 

Jo 1 + $(t) 2 



exp(-2t 2 + tx 2 2v / 2) ^ exp(-t 2 + tx 2 2y/2) 

\- y'Ttt — 



-An 3 2V2vr / dt 



-A 2 n 4 3V2V 



df 



3 

\ 2 2 v - 

—A n — exp 
6 



,y/2 



dtt 



-2x1) 



exp(-3t 2 + tx 2 2V2) 

[i + w + 

exp(-4t 2 + ta; 2 2v / 2) 

[1 + W 
exp(-3t 2 + te 2 2V2) 



vrt 



1 + 



exp(-2t 2 + te 2 2V2) 



+ 0rt 



[1 + «&(*)]* 



+ V^t 



[1 + $(t)] 2 

exp(-3t 2 + tx 2 2^2) 
exp(-2t 2 + tx 2 2y/2) 



[1 + *(*)]= 



+ -A 2 n 



24 9a; 2 



x < exp(— 2x 



2 

oo 



dt 



°° exp(-2t 2 + tx 2 2y^) 
[l + $(t)] 2 

X2 <9x 2 / 
exp(-2t 2 + tx 2 2y/2) 



2 ^ ^ -n 2 v / 2exp(-2a; 2 /A 2 ^- — ^ 



exp(-t 2 +to: 2 2v / 2) 
1 + $(t) 



[1 + $(t)p 



+ 



1.2) 



On the other hand, n qu (x 2 ) is given by (I6.2p . 

The classical part of the sum rule comes from the second line of (18.21) and the first 
line of (|B2|). Since exp(-t 2 )/[l + $(t)] 2 = -( > /7r/2)(d/d*)[l + $(£)] _1 , an integration 
per partes gives, with n = l/n, 

1 



y dxiXi J dt/in (2)T (ri, r 2 ) + x 2 n(x 2 ) 



exp(— 2x 2 ). 



•3) 



v^vnV 2 

Integrating on x 2 gives the classical result — l/(47r). 

Using (I6.3p and (18. ip . it is straightforward to integrate on x 2 most of the terms of 
(18. 2p and x 2 times (16. 2p . The term in (18. 2p . line —3, of the form d 2 F(x 2 )/dx 2 requires 
some care for evaluating its integral on x 2 , which is dF(x 2 )/dx 2 |g° . The point is that 
F(x 2 ) has a term which is not zero at infinity. Indeed 

exp[-(t-W2) 2 ] 



exp (—2x2 



°° d ^exp(-t 2 + tx 2 2V2) 



o 



dtt- 



1 + 
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>V2 r dt' eM ; f2) = ^/2 X2 (8.4) 



~ £ 2 1 



where t' = t — x 2 v / 2, and the derivative with respect to x 2 of ( 18.4ft is y 7r/2 at infinity. 
The contribution of this derivative at infinity to the sum rule is found to be — A 2 /(247r). 
Similarly, the term in ( 16.21) . line —2, gives a contribution of the form X2[d 2 G(x2)/dx 
when integrated per partes on x 2 , in particular it gives a term —G(oo). Since 



1]; 



V ' Jo l + $(t) x^™J-oc 2 V 7 ' V ^ 

the corresponding contribution to the sum rule is — A 2 /(247r). All other quantum 
contributions to the sum rule are found to cancel each other, at order A 2 . Finally, 

-i-&« +0 ^- (8 - 6) 

Since A is in units of a, the sum rule ( 14.51) does have a quantum correction 
A 2 /(127ra 2 ) at order A 2 . 



oo 

dxo 



r 

dxtxt / d?/i^ 2)T (ri,r 2 ) +a; 2 n qu (x 2 ) 



9. Conclusion 

That the perfect screening sum rule (14. ip is satisfied is no surprise. This is rather 
a check of our calculations. The other sum rules are less straightforward. In their 
heuristic macroscopic derivations [8], an essential feature of the quantum systems, 
that the n-body densities have to vanish on the wall, was not explicitly taken into 
account. Thus a check that they are indeed correct is welcome. 

The case of the asymptotic-form sum rule (14.31) is not entirely satisfactory. Although 
we have checked it at order A 2 , at this order we could only verify that there is no quantum 
correction, in agreement with the expansion of the rhs of (14.31) . It would have been more 
satisfactory to verify the quantum correction of order A 4 of this rhs. Unfortunately, this 
would involve pushing the Wigner-Kirkwood expansion to that order, which would be 
straightforward but so tedious that we cannot hope for its feasibility in a reasonable 
time. 

The dipole sum rule (14. 5 p is more tractable, and we have indeed checked a finite 
quantum correction at order A 2 . Moreover, another derivation of this sum rule is feasible. 
Indeed, the generalization of the classical Stillinger-Lovett sum rule [TO] for a bulk 
quantum one-component plasma has been microscopically done [UJ; adapting the result 
to two dimensions gives, with r = r x — r 2 , 

/ drr^'L^r,,,,) = --L^coth^p. (9.1) 

Then, we can use the same kind of argument as the one in pages 965 and 966 of ref. [8] 
for deriving the dipole sum rule (14.51) . 
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